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Abstract. We discuss a detailed phase diagram and other microscopic characteristics 
on the applied magnetic field - temperature (H a — T) plane for a simple model of 
correlated fluid represented by a two-dimensional (2D) gas of heavy quasiparticles 
with masses dependent on the spin direction and the effective field generated by 
the electron correlations. The consecutive transitions between the Bardeen-Cooper- 
Schrieffer (BCS) and the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phases are either 
continuous or discontinuous, depending on the values of H a and T. In the latter case, 
weak metamagnetic transitions occur at the BCS-FFLO boundary. We single out two 
different FFLO phases, as well as a reentrant behaviour of one of them at high fields. 
The results are compared with those for ordinary Landau quasiparticles in order to 
demonstrate the robustness of the FFLO states against the BCS state for the case with 
spin-dependent masses (SDM). We believe that the mechanism of FFLO stabilization 
by SDM is generic: other high-field low-temperature (HFLT) superconducting phases 
benefit from SDM as well. 
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1. Introduction 

Unconventional superconductivity in heavy-fermion and organic-metal systems is 
studied almost as frequently as high-temperature superconductivity and comprises a 
number of heavy-fermion and organic metallic systems pQ. Among the states observed 
and discussed intensively recently is the superconductivity in the systems without space 
[2], and time [3] inversion symmetry, the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state 
[I] and the states in which magnetic order, usually antiferromagnetic (AF), coexists 
with the FFLO [3 [6] or the Bardeen-Cooper-Schrieffer (BCS) type of state. Also, 
the FFLO states are discussed recently in the context of cold atomic fermionic gases 
[7] and quark-gluon plasma [8]. One of the basic motivations for these studies is the 
circumstance that the normal state can be represented by a Fermi fluid, albeit almost 
localized, so the nature of paired state can be rationalized to a greater detail. Also, the 
intriguing feature of those superconductors is a cooperation rather than competition 
with magnetism [9]. The field- induced magnetism can be incorporated into the FFLO 
state, since there are substantial portions of the Brillouin-zone volume over which the 
quasiparticle excitations are gapless. 

It is the later topic (the FFLO appearance) which is the principal subject of 
this paper starting from a two-dimensional (2D) d-wave superconductor composed of 
unconventional (correlated) quasiparticles. Namely, we represent the heavy-fermion 
liquid by a gas of quasiparticles with the spin-direction dependent effective masses 
(SDM), which were indeed observed in CeCoIns and other systems pU] (in that case 
m* = m a , o = ±1 being the particle spin quantum number). Another non-trivial feature 
of our approach is the inclusion of the effective field h cor acting upon the magnetic 
moments in the spin-polarized state [IT] . Both characteristics are generated by the 
electron correlations treated in the mean- field-type schemes [TT| [12] . The experimental 
motivation for our study is the observation of both SDM \X0\ and FFLO (or FFLO mixed 
with magnetism) [5] in the same heavy fermion system CeCoIn 5 , and the question we 
tackle is whether these two phenomena are interconnected. To address it we consider a 
simplest situation of electron gas with the FFLO state in the simplest form (FF type 
with A(r) = A e iqr ). We also consider a c?-wave form of the superconducting gap, 
i.e. Ak,Q = Ac^cosfca; — cosk y ), where Q is the Cooper pair momentum (Q / 
in FFLO state). Such form of the gap reflects the principal feature of quasi- two- 
dimensional superconductivity in strongly correlated electrons [15] . We show that the 
phase diagram with these high-field low-temperature (HFLT) c?-wave superconducting 
phases in 2D differs remarkably from its 3D correspondant with the s-wave symmetry 
[TH H3) . Namely, several FFLO states appear in the present situation even when we 
disregard the possibility of their coexistence with antiferromagnetism (5J [6]. We also 
show that a weak metamagnetic transition accompanies the BCS — > FFLO discontinuous 
transition. 
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2. Model: unconventional gas of quasiparticles with real-space pairing 

The principal features of our approach have been defined earlier [H] (cf. particularly 
Sections V - VII). We start from the effective quasiparticle picture which is common 
to both narrow-band and hybridized correlated-electron systems [HI [El El [17] . The 
explicit form of quasiparticle energies in the gas of correlated quasiparticles is 

h 2 k 2 

£k<x = — 0-(h + hear) - fl, (1) 

where /x is the chemical potential, h = g[iBH a , and the mass enhancement factor 
in the large-t/ limit [17] in the simplest situation is 



m a 1 — n CT l — n/2 m _ 

m>B 1 — n 1 — n 2(1 — n) 

= — (m av -aAm/2), (2) 
m B 

with fn = n-j- — n± being the spin polarization and n = n-f- + the band filling. 
Also a — ±1 is the spin quantum number, is the band mass, Am = ni2 — m\ is the 
mass difference and m av = (mi + rri2)/2 is the average mass. As one can see, the spin- 
dependent mass enhancement is particularly strong for an almost half-filled case when 
1 — n = 8 1, i.e. for the quasiparticles close to the Mott-Hubbard localization. Here 
the superconducting phases in this 2D <i-wave superconductor are discussed in detail 
and compared briefly with the previous results [HI [15]. In connection with this one 
should note that the concept of SDM has been also used in the context of coexistence 
of ferromagnetism and superconductivity [18] . 

Even though our considerations represent a model situation, we assume the 
following values of the parameters, emulating the heavy fermion systems: the filling 
n = 0.97, the elementary square-cell area S = (4.62A) 2 , the starting (H a = 0) 
quasiparticle mass m av = 100m (data for CeCoIn 5 [TO]), the pairing potential cutoff and 
magnitude hue = 17 K, and Vo = 90 K, respectively. The characteristic energy scale 
associated with spin-fluctuations in CeCoIn 5 is T s f = 10 K [19] - a value comparable to 
our fuvc- F° r those parameters, the chemical potential was equal to \i ~ 126 K. This 
means that Vq < \i and the (weak-coupling) BCS approximation can be regarded only 
as a proper solution on a quantitative level at best. Additionally, the chemical potential 
is readjusted in the superconducting state so that n is constant. The pairing potential 
has the separable <i-wave form 

Vk,k' = — Vo(cos k x — cos hy)(cos k' x — cos k' ), (3) 

which differs slightly from that used in Ref. [l_5]d. For thus defined quasiparticles 
with energy £ kcr , we derive their correspondants E ka in the superconducting states 

[HIED! 
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E ka = E k + a&\ (4) 
£ k = \A£ )2 + A£ iQ , (5) 
£k = 2^ k +Q/ 2 t + ^-k+Q/24), (6) 

Ck""* = 2^ k +Q/ 2 t ~~ C-k+Q/24.), (7) 

as well as the free-energy functional J 7 and the system of four self-consistent integral 
equations for the field h cor , magnetization m, gap magnitude Aq and the chemical 
potential //. Explicitly, starting from the free energy functional J 7 , we obtain the 
corresponding integral equations of the following form 



7 = -^T^ln(l + e-^) + ^(^ ) -E k ) + 



kcr 

.2 



An iV 

+ N-^- + /iN+—mh cor , (8) 
Jy ^— ' am 

kcr 



??? 



N ^ dm \ E k 
k K 

kcr 



Aq = ^(COS ^ - COS 



2 X 

k 



1 - /(£ kt ) - f(Ek^) 



Aq, (113 



2£ k 

n = n-j- + = 

= ^ EK/^) + ^ - ( 12 ) 

kcr 

Those quantities determine the physical free energy in different (BCS, FFLO, 
NS) states which are compared to obtain the phase diagram and other microscopic 
characteristics, as we discuss below. Note that we limit ourselves to a single Q (Fulde- 
Ferrell type) solution [HI [20] as we intend to describe superconductivity with SDM in 
the simplest case and thus test the importance of the quasiparticles mass spin-direction 
dependence. In that situation, the whole problem comprises a simultaneous solution of 
those four integral equations for fi, m, Aq and h cor for fixed Q followed by subsequent 
minimization of thus obtained physical free energy T with respect to Q. 
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Figure 1. (Colour online). Phase diagram for the cases with the spin-dependent 
(SDM), a), b) and with spin-independent masses (SIM), c). Light (yellow) regions 
correspond to Q = (BCS phase), the darker (blue, red) to the state with Q ^ 
(FFLO phase) and the white to normal state (NS). The colour scale in a) is defined 
by the pair momentum Q. Note that for SDM with increasing temperature, the 
transition from BCS to FFLO state occurs at higher fields, in qualitative agreement 
with experimental results [5]. The different FFLO phases are exhibited in b). The 
red region corresponds to the Cooper- pair momentum Q in the k x direction (0q = 0), 
whereas the blue one to the momentum along the diagonal (k x = k y , 9q = 7r/4). Note 
that this anisotropy results solely from the d-w&ve gap symmetry, as the unpaired gas 
is isotropic. The dashed line marks the BCS critical field H C 2 in the Pauli limit |21) . 
and the dot-dashed line marks H C 2 for the solution with 6q = 0. 
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Figure 2. (Colour online), a) Gap parameter Aq (in units of K) and b) Cooper 
pair momentum Q in units of the Fermi momentum difference Akp = fep T — /cfj,, 
both on H a - T plane. Transitions between various phases are seen as a change of the 
magnitude of the gap: the lower-field transition are first-order, whereas the transition 
to normal state is continuous. 



3. Results: BCS vs FFLO states 

The overall phase diagram in 2D case on the applied magnetic field (H a ) - temperature 
(T) plane is exhibited in figure [H for the cases with spin- dependent (SDM) (a, b) and 
the spin-independent (SIM) (c) effective masses. The FFLO phase is robust only in the 
former for the s-wave solution for the three-dimensional gas [14], although the 

difference is greater in the 3D case. The specific difference is that in the present case 
two distinct phase-boundary lines appear inside the FFLO state, as detailed in Figure 
lb: the topmost and the lowest parts (red colour) have the Cooper-pair momentum 
Q oriented along the k x (or k y ) direction, whereas the middle phase (blue colour) has 
Q along the diagonal ik x = k y ). Also, superconductivity of FFLO type exists up to 
the field of 35 T in the SDM case, i.e. the field more than 4 times larger than that 
for the SIM case. Hence, the former system indeed belongs to the class of high-field 
low-temperature superconductors. 

To visualize the detailed nature of the transition to the FFLO phase we have plotted 
in figure [2] profiles of the gap magnitude Aq and the Cooper pair momentum |Q|, both 
on the H a - T plane. In the low-T limit the observed gap jumps meaning that the 
transitions BCS — > FFLOl (Q || k x &xis), as well as the transition FFLOl — > FFL02 
(Q || (k x , k y ) diagonal) and FFL02 — > FFLOl' (Q || k x axis) are discontinuous, whereas 
the transition to the normal state is continuous (cf. also [12J). As the temperature 
increases, all the transitions (except that from FFL02 to FFLOl') become continuous, 
but the exact position of the terminal bicritical point will not be discussed in detail 
here. The phase FFLOl' illustrates a reentrant high-field behaviour for FFLOl phase 
(cf. also figure [lb for the SIM case). Note also that the FFLO states exist far beyond 
the second critical field H c2 |21j for BCS state, marked by the dashed line. 

The above phase transitions can be connected with the magnetization changes. 
This is because the FFLO phase encompasses semimacroscopic regions of k-space with 
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Figure 3. (Colour online). Spin polarization m = n-\—n\. as a function of applied field. 
Note the weak jumps corresponding to the discontinuous transitions at T = 0.02 K 
and T = 0.50 K for SDM case (a) and much larger in the SIM case. (b). For the SDM 
case all transitions at T = 1 K are continuous. 



gapless quasiparticle excitations in the superconducting phase. This means that the 
magnetization curve will show a nontypical behaviour, particularly in the vicinity of 
the transition to FFLO state, as displayed in figure El Namely, the m(H a ) exhibits 
a weak metamagnetic behaviour accompanied by a weak jump at the two lower-field 
transition points. It is surprising at first look that the corresponding jump is much 
larger in the SIM case. However, one must remember that in the SDM case the field h cor 
compensates largely the applied field (see figure 0] for details). The spin magnetization 
does not include the magnetic dipole moment which may arise from an inhomogeneous 
current-carrying state when Q ^ 0. Obviously, the FFLO state may coexist with AF 
(or SDW) or spin-flop phases, but these cases are not discussed here, as we would like 
to characterize in detail here the "pure" FFLO state to single out its novel features in 
the SDM case. 

To compare our results with those for three-dimensional system and s-wave pairing 
symmetry we recall here the mechanism behind the FFLO stabilization by SDM 
presented in [H] (cf. Section VI there). Namely, SDM compensate the Zeeman effect 
influence by reducing the Fermi wave vectors splitting. Therefore, superconducting 
state with SDM has higher critical fields (here h C 2 = 10 T for SIM, and h C 2 = 36 T 
for the SDM case, cf. figure [T|). The FFLO state benefits from SDM by a greater 
extent than BCS because spin polarization m in the latter is smaller (cf. figure |3]), 
and from (J2J) the mass difference Am oc m. Therefore, in BCS the mass difference 
is smaller, and the Fermi wave vectors splitting larger than in FFLO (the Zeeman 
term influence is compensated less effectively). Hence, at T = the FFLO fills about 
1/2 of the phase diagram for SIM, and about 2/3 for SDM. On the other hand, as 
temperature T increases, the spin polarization increases in the BCS state (see figure [3]) 
allowing larger mass difference Am and reducing Fermi wave vectors splitting enhancing 
superconductivity. This is why the transition line between BCS and FFLO is curved 
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Figure 4. (Colour online). Correlation induced effective field h cor in the 
superconducting state on H a — T plane relative to that in the normal state h cor FSi 
which is typically equal to — 0.33/i, and therefore the field h cor compensates applied 
field H a . 

upwards in the SDM case. In the present situation, the BCS state can have a substantial 
spin-polarization already at T = (unlike in the 3D, s-wave case) and therefore the BCS 
state can benefit from SDM already at T = 0, and the FFLO state is not stabilized so 
spectacularly here, as it was in the 3D case (where in the BCS phase fn ~ at T = 0). 

For the sake of completeness, we draw in figure H] the effective field induced by 
the correlations. The jumps reflect the discontinuous transitions discussed above. The 
field h cor (in units of h cor ps f° r the unpaired Fermi sea) increases both with increasing 
temperature and field. The mass difference Am = — m-j- changes with the applied 
field reflecting the change in fn(H a ); the relative difference Am/m av reaches about 10% 
for the applied field of the order of 30 T. 

4. Conclusions 

In summary, we have singled out different FFLO states in a 2D gas of correlated 
quasiparticles with spin-dependent effective masses (SDM) and effective field induced by 
the electron correlations, as well as compared them briefly with those in SIM case [15] . A 
number of FFLO phases appears and these phases are stable in an unusually high fields 
only for the case with SDM which were indeed discovered in CeCoIn 5 and other systems 
[10] . It is suggested that these nonstandard properties of quasiparticle states should be 
their universal feature for all the systems close to the /- or <i-electron Mott-Hubbard 
localization if the atomic disorder effects are very weak. Namely, other HFLT phases 
(including various FFLO phases mixed with antiferromagnetism) can be stabilized from 
having SDM as well, since they always have higher spin-polarization than the uniform 
superconducting state, and then SDM compensate the Zeeman term influence more 
effectively than in the uniform superconducting state. Extension of these results to 
incorporate the antiferromagnetic ordering within the present approach and for realistic 
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highly anisotropic three dimensional electronic structure, would most probably provide 
a decisive answer about the nature of high-field low-temperature phase in CeCoIn 5 [5] 
and organic systems [22] . 
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